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We investigate the effect on high-order harmonic generation of the distortion of molecular orbitals 
by the driving laser field. Calculations for high-order harmonic generation including orbital distor- 
tion are performed for Na. Our results allow us to suggest that field-distortion is the reason why 
the two-center interference minimum has never been observed experimentally in Na. We propose 
experimental parameters which should allow an observation of the two-center interference minimum. 

PACS numbers: 33.20.Xx,42.65.Ky,42.50.Hz 



I. INTRODUCTION 



During the last decades, high-order harmonic genera- 
tion (HHG) has been subject to extensive experimental 
and theoretical studies. Apart from being a source of co- 
herent XUV radiation [TJ [5] , and ultrafast pulses [3HS] , 
the process of HHG has also been used to extract in- 
ternuclear separations [HHS] and electronic orbitals jTDJ- 
[2]. The intrinsic timing in the HHG process associ- 
ated with the time from ionization to recombination can, 
e.g., be used to resolve nuclear dynamics on an ultra- 
fast timescale [6 8J. Turning to orbital tomography, the 
HHG signal has been used to reconstruct molecular or- 
bitals: the highest occupied molecular orbital (HOMO) 
in N 2 [3I5HI3!, the HOMO-1 in N 2 13 , the HOMO in 
C0 2 [H 03] and the HOMO in 2 \l2\. Given the ex- 
perimental laser intensities, it is surprising that, in all 
of these cases, there has been no measurable signature of 
field-distortion in the reconstructed orbitals. The present 
work identifies for the first time a clear signature of sig- 
nificant field-induced orbital distortion in the HHG pro- 
cess by focusing on the effect of the distortion on the 
two-center interference minimum [15J-T8] . Indeed, one 
of the most prominent features of an HHG spectrum is 
the presence of a minimum. Such a minimum may have 
different origin depending on the target and laser param- 
eters, but it always carries important information about 
the target. For example, a Cooper minimum holds in- 
formation about the electronic structure of the target 
[I9H23] . Participation of multiple orbitals in the HHG 
process gives rise to a so-called dynamic minimum which 
allows elucidation of the interplay between the HOMO 
and lower-lying orbitals in the HHG process [2H US] . Fi- 
nally, a two-center interference minimum gives access to 
structural information [TotilSj . 

Accounting for the effects of orbital distortion allows 
us to resolve a long-disputed puzzle in N 2 : Both exper- 
imental and theoretical works have proven the existence 
of Cooper minima [T9H23] [26] and the role of multiple 
molecular orbital contributions in HHG spectra [57] . Yet 
the two-center interference minimum, though theoreti- 
cally predicted [28j [29] , has never been observed experi- 
mentally [10l EH [30] • We explain the lack of a two-center 
interference minimum in terms of laser-induced orbital 



distortion. We also show that for specific field parame- 
ters and orientation of the molecule it should be possible 
to observe experimentally the two-center minimum in the 
long trajectories of the HHG spectrum. 



II. MODEL 

We assume that the electronic structure of the 
molecule adiabatically adjusts to the time- varying exter- 
nal field, F(t). With this assumption the Lewenstein 
model [3.1] is extended to include not only the Stark 
shift [32], but also the field-distortion of the initial and 
final states. For sake of clarity we present the approach 
with participation of only the HOMO in the HHG pro- 
cess (the methodology can be straightforwardly extended 
to include more orbitals). The HHG spectrum along the 
polarization e can then be calculated as |33| 



S e (oj) = 



(«(*)) 



z iu;t dt 



(1) 



where the expectation value of the dipole velocity oper- 
ator reads 



(«(*)> = ~i 
with 



«!/' / AkM, cc {t)e- %s ^^ M ion {t')+c.c., 

(2) 



M ion (t') = (^(t')\V L (m(F(t'))) 
the ionization matrix element, 



(3) 



S(k, t, t') = J l - (k + A(t")f dt" - J E 



{t")dt" (4) 



the phase accumulated while the electron propagates 
from time t' to t, and 



M iec (t) = (ip(F(t))\v\iljl(t)), 



(5) 



the recombination matrix element. The Stark-shifted 
HOMO energy is denoted by E(t) in Eq. Q. Using per- 
turbation theory E{t) = E -fJ,- F(t) - l/2F(t) T g±F(t) 
with Eq the field-free HOMO energy, fi the dipole of the 
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HOMO and a the polarizability tensor. In Eqs. ^ and 
([5]) the HOMO \ip(F(t))) depends explicitly on the in- 
stantaneous value of the driving laser field F(t). Propa- 
gation of the continuum electron with momentum k is de- 
scribed by the Volkov state \ip~l (*)), and V L (t) = F(t) ■ r 
describes the interaction with the laser field. We ob- 
tain initial and final states by calculating field-dependent 
orbitals with quantum chemistry methods [34] using an 
aug-cc-pTZV basis set [35] |3J>] • This approach was suc- 
cessfully used to study deformations in CO2 molecules 
subject to even stronger fields than used here [37]. To 
obtain the spectrum, the field-distorted orbitals and cor- 
responding energies are calculated at every time step 
needed for an accurate evaluation of the integral in 
Eq. ([2]). This approach incorporates many-electron ef- 
fects on the bound state since the calculation of molec- 
ular orbitals exposed to the laser field involves the mul- 
tielectron response. We note that the present approach 
requires that the distorted orbitals are calculated at the 
instantaneous value of the external field. Hence, the ap- 
proach refers to a description using the length gauge 
for the electron-laser interaction, which allows the us- 
age of standard quantum chemistry program packages 
for static-field calculations. The HOMO obtained in this 
way may be transferred to the velocity gauge by a suit- 
able unitary transformation [55] , 



III. RESULTS 

We freeze the nuclei in N2 at the internuclear distance 
of 2.08 a.u. and set the ionization potential to 15.56 eV. 
The molecule is aligned [39] at j3 = 75° with respect to 
the polarization axis of a linearly polarized driving pulse. 
The pulse has a five-cycle trapezoidal envelope with one 
optical cycle for ramp-up and -down. The peak inten- 
sity is 9 x 10 13 W/cm 2 , and the carrier wave length is 
1300 nm. For this orientation and intensity it can be 
safely assumed that multiple orbital contributions are 
negligible (see Ref. [27] ) and therefore we restrict our- 
selves to the signal emitted by electrons detached from 
the HOMO. Possible depletion effects (not treated here) 
are not expected to affect the shape of the HHG spectra, 
but only the overall yield, and will, hence, not alter our 
main conclusions. 

When field-free HOMOs are used in the ionization 
[Eq. ([3])] and recombination [Eq. |5])] steps the model 
reduces to that of Ref. [32], i.e., the Lewenstein model 
[3T] including Stark shifts. This model predicts a two- 
center interference minimum at - 58 eV [Fig. Qa)]. The 
minimum is identified as a two-center interference min- 
imum by its change in position when varying the align- 
ment angle: We have checked that the minimum moves to 
higher photon energy when increasing the alignment an- 
gle [15l[l6l [29]. Such a minimum has never been observed 
experimentally (see, e.g., Refs. [TO] [30]), irrespective of 
orientation and laser intensity. The two-center interfer- 
ence minimum of concern should not be confused with 
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FIG. 1: (Color online) HHG spectra for N2 (see text for laser 
parameters and orientation). The thick (red) curves indi- 
cate spectra smoothed with a Gaussian distribution with full- 
width-half-maximum of 4 harmonics, (a) Field-free HOMO 
in ionization [see Eq. |3|] and recombination [see Eq. |5]l] 
steps, (b) Field-distorted HOMO in ionization and field-free 
HOMO in recombination steps, (c) Field-free HOMO in ion- 
ization and field-distorted HOMO in recombination steps, (d) 
Field-distorted HOMO in ionization and recombination steps. 



the Cooper minimum, which appears at ~ 38 eV [53J [2S] • 
The latter minimum is related to scattering properties of 
the molecular potential [40] and is not considered in the 
present work, where we focus at larger energies, where the 
approach treating the continuum states as Volkov waves 
is reasonably accurate. Next we use a field-distorted 
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HOMO in the ionization matrix element, leaving the 
state in the recombination matrix element field-free. The 
two-center interference minimum remains visible in the 
spectrum [Fig. [ljb)], and hence we reach the conclusion 
that the deformation of the initial state does not signifi- 
cantly obscure the structural information associated with 
the two-center minimum. An important change as com- 
pared with Fig.[]Ja) is a much lower yield which is caused 
by destructive interference between amplitudes originat- 
ing from different ionization times. In the next step we 
examine the role of the recombination matrix element [see 
Eq.Q] and use a field-distorted HOMO in the recombi- 
nation, and a field-free HOMO in the ionization step. 
The result is depicted in Fig. [TJc) . The disappearance of 
the minimum proves that the recombination matrix ele- 
ment is sensitive to orbital distortion, and in this sense 
holds structural information about the target. Similarly, 
the minimum disappears when taking into account the 
field-distortion of the HOMO in both the ionization and 
recombination steps [Fig.JlJd)]. 

The two-center interference minimum appears when 
the recombination matrix elements on both atomic cen- 
ters are of similar magnitude and their phase difference 
is 7r. To calculate the recombination matrix element on 
a given center we use that the orbitals are expanded in 
atomic centered Gaussian basis functions and take the 
part of the molecular orbital centered on the atomic cen- 
ter of interest |41| . In Fig. [2] we show how orbital dis- 
tortion significantly changes the relative strengths and 
phases. If the HOMO is field-free in the recombination 
step [Figs. [lja) and (b)], then the norms of recombina- 
tion matrix elements on both centers are equal as shown 
in Fig. [2ja). The phase difference reaches tt for a con- 
tinuum electron energy of ~ 42 eV corresponding, by 
adding the ionization potential, to a harmonic energy 
of ~ 58 eV. This explains the position of the minima 
in Figs, [ija) and (b). When taking into account field- 
distortion in the recombination step the value of the field 
changes with recombination time. In Figs. 2Vb) and (c) 
we depict the norm and phase of the recombination ma- 
trix elements when the HOMO is distorted by F = 0.01 
a.u. and F — 0.04 a.u., respectively. These figures show 
that the norms on different centers become increasingly 
different for increasing field. As a consequence the two- 
center interference minimum weakens and subsequently 
totally vanishes. Monitoring the variation of the mini- 
mum with the field strength allows us to specify a range 
of field values for which observation of the minimum is 
possible. This range corresponds to recombination at 
fields \F\ < 0.012 a.u., and we will refer to it as the 
weak-field region. 

While Fig. [2] defines the weak-field region, it does 
not contain information about the value of the field at 
the time of recombination. To this end we consider a 
time-profile analysis [42H44] . Figure [3] presents such an 
analysis corresponding to the spectra in Figs. 1(a) and 
1(d). The figure shows which recombination, i.e., emis- 
sion times contribute to the signal at a given photon en- 
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FIG. 2: (Color online) Norm and phase (inserts) of recom- 
bination matrix elements [Eq. |5]l] for (a) F = a.u., (b) 
F = 0.01 a.u., (c) F = 0.04 a.u. The full (blue) curves show 
results using the field-distorted HOMO, dot-dashed (red) 
curves denote field-distorted recombination on the N atom 
at the negative z axis in the molecular frame, and dashed 
(green) denote field-distorted recombination on the N atom 
at the positive z axis in the molecular frame. 



ergy. Figures 3(a) and (c) show that the harmonics come 
in bursts each half cycle. Each harmonic burst consists 
of two branches, known as the short and long trajecto- 
ries [45]. The short trajectories start less than 2/3 of 
an optical cycle before recombination and give rise to 
branches with positive slope. The long trajectories start 
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FIG. 3: (Color online) Emission times of the spectra presented in Figs. 1(a) and 1(d). (a) Emission times when performing 
the calculation with undistorted HOMOs in the ionization and recombination steps [Fig. 1(a)]. Curves indicate positive (full, 
white) and negative (dashed, white) values of the laser field at the recombination time, (b) Field-free, undistorted HOMO, 
(c) Emission times when performing the calculation with field-distorted HOMOs in the ionization and recombination steps 
[Fig. 1(d)]. Curves indicate positive (full, white) and negative (dashed, white) values of the laser field at the recombination 
time. The vertical (yellow) lines identify the weak-field regions (see text), (d) Field-distorted HOMO for F = 0.04 a.u.. (e) 
Field-distorted HOMO for F = 0.01 a.u. In (d) and (e) full (black) arrows indicate the field direction. The arrows from (d), 
and (e) to (c) indicate the recombination times for which the HOMOs are depicted, and dashed (black) lines in (c) guide the 
eye to the corresponding value of the field. 



more than 2/3 of an optical cycle before recombination, 
and give rise to branches with negative slope. Figure 
3(a) shows that when the calculation is performed us- 
ing the undistorted HOMO shown in Fig. 3(b), there is 
a minimum at around 58 eV in both branches, consis- 
tent with the minimum in Fig. 1(a). Figure 3(c) shows 
that when field-distortion is included (see Figs. 3(d) and 
3(e) for snapshots of the distorted orbitals), the short 
trajectory branches are more intense than the long tra- 
jectory branches, and no minima occur in the former. For 
example, the short trajectory branches give rise to pho- 
tons at 61 eV while the corresponding long trajectory 
branches do not contribute significantly to this photon 
energy. In Fig. 3(c) the vertical (yellow) lines indicate the 
time windows where the field at recombination is within 
the weak-field region. We see that the trajectories that 
rccombinc in these windows belong to the long trajec- 
tory branches. The minimum at 61 eV for the long tra- 
jectory branch is the two-center interference minimum, 
that is shifted in position by orbital-distortion from 58 
eV in the undistorted case [Fig. 3(a)]. Following the dis- 
cussion in connection with Fig. 2, this implies that the 
two-center interference minimum can only be observed in 
the long trajectory branch, and therefore a measurement 
isolating the contribution to the spectrum from the long 
trajectories may lead to an observation of a two-center 
interference minimum in N 2 . 

The reason why the minimum is not seen in Fig. [l|d) 
is that the short trajectory harmonics dominate the spec- 



trum in our calculation (see Fig. [3]). To isolate the spec- 
trum originating from the long trajectories, macroscopic 
phase matching should be used to select the long trajec- 
tories by either changing the focusing conditions of the 
driving field [35] or by far-field separation jJTj (see also 
the review [45] and references therein). 

IV. CONCLUSION AND OUTLOOK 

Our results indicate that field-distortion of molecular 
orbitals can be important for understanding spectral fea- 
tures of molecular targets. The effect of orbital distor- 
tion is manifested in both ionization and recombination 
events. The orbital distortion in the ionization event 
results in a lower HHG yield in N 2 compared to the 
undistorted case. Orbital distortion in the recombina- 
tion step influences the structural information encoded 
in the spectrum. For example, the absence of a two- 
center interference minimum in the HHG spectrum of 
N 2 is a clear signature of orbital distortion taking place 
when the short trajectory recombines outside the weak- 
field limit. By carefully adjusting the molecular orien- 
tation and the laser parameters, conditions can be ob- 
tained such that the long trajectory recombines within 
the weak- field limit, allowing the observation of the N 2 
two-center interference minimum. 

The field-distorted orbitals can be obtained in stan- 
dard quantum chemistry software packages [34], and 
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therefore the methodology can be applied to many molec- 
ular systems. The approach opens up the possibility to 
study a range of questions that can shed more light on the 
process of HHG. For example, the effect of the orbital dis- 
tortion will be very sensitive to the value of the polariz- 
ability. N2 has a relatively high polarizability, and there- 
fore the effect on the spectrum is large. In other systems, 
for example CO2, the polarizability is much smaller, and 
orbital distortion should be less important. In addition, 
an approach including orbital distortion is highly rele- 
vant for orbital tomography. Our work indicates that 
HHG-based orbital tomography will never allow the re- 
construction of a field-free orbital, but rather gives a co- 
herent average over field-distorted orbitals. The degree 



of field-distortion may be reduced by using short driving 
pulses [48], and isolating the contribution from the long 
trajectories. Finally, we note that the present approach 
can also be used to address the effect of orbital-distortion 
in the strong-field approximation for ionization [49H51) . 
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